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We show that when non-commutative quantum mechanics is formulated on the Hilbert space of 
Hilbert-Schmidt operators acting on a classical configuration space, spectral triplets as introduced 
by Connes in the context of non-commutative geometry arise naturally. A distance function as 
defined by Connes can therefore also be introduced. We proceed to give a simple algorithm to 
compute this function in generic situations. Using this we compute the distance between pure and 
mixed states on quantum Hilbert space and demonstrate a tantalizing link between statistics and 
geometry. 
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I. INTRODUCTION 

The idea that space-time may be quantized in some sense and that a more abstract notion of 
geometry may be required has become firmly established in the physics literature in the past two 
decades. Connes was the first to introduce, in a rigorous fashion, the notion of non-commutative 
geometry [![, captured by a spectral triplet (A, H, D) with A an involutive algebra acting on a Hilbert 
space T~L by a representation 7r and D is the so-called Dirac operator on H. The further technical 
conditions that this operator must satisfy can be found in 

In a separate, yet related development, it was argued by Doplicher et al [2] from the considerations 
of both general relativity and quantum mechanics that the localization of an event in space-time 
with arbitrary accuracy is operationally impossible and this feature is captured by postulating a non- 
vanishing commutation relations between operator-valued coordinates. In its simplest form they are 
given as 

where is anti-symmetric and its entries are viewed as new fundamental constants. This form of 
non-commutativity also follows from the low energy limit of string theory Q. 

Despite the fact that the two dimensional non-commutative (Moyal) plane also fits in the framework 
of non-commutative geometry as was established quite some time back for both compact Q and non 
compact cases [5|, a detailed and explicit investigation into the link between these two versions of 
non-commutative space-time and associated geometries seems to have been started only recently with 
explicit computations of the Connes' spectral distance between pure state harmonic oscillator functions 
as well as translated states [f| 0] • One possible reason could be that the explicit computation of the 
Connes' distance function can be very challenging. Another reason may be that the Connes' approach 
was primarily applied to " almost commutative spaces " which are generically of the form M x F with 
M representing 4— dimensional ordinary commutative spacetime describing the gravitational part and 
F corresponding to a zero— dimensional discrete space described by matrices, which takes care of 
the gauge part of the standard model [1]. The spaces of the type M x F has thus only a " mild " 
non-commutativity stemming from this internal space F. This is believed to be adequate to describe 
physics up to the GUT scale. However, at a still higher energy scale, like in the vicinity of the Planck 
scale, one may have to take the fuzziness of spacetime, captured by equations like |T]), seriously. 

Connes' spectral distance function relies on the notion of states, which are positive linear functionals 
of norm one, over the involutive algebra A. They are also closely related to the notion of density 
matrices, which describe the physical states of a quantum system. This naturally raises the question 
whether the formulation of non-commutative quantum mechanics as set out in fiol . [IlT | , which entails 
a representation of the quantum system on the Hilbert space of Hilbert-Schmidt operators, does not 
naturally encode geometric information in the context of Connes' non-commutative geometry. This is 
the question we investigate here and we indeed find that this formulation of non-commutative quantum 
mechanics naturally encodes geometry in the form of spectral triplets. Moreover we succeed in giving 
a simple algorithm to compute the Connes' spectral distance between states. It should, however, be 
remarked that the notion of spectral triplets that we introduce here is stronger than the conventional 
one, which takes the involutive algebra to be the space of bounded operators on %, while here it 
is taken as the Hilbert-Schmidt operators and therefore is itself a Hilbert space. It is precisely this 
Hilbert space structure of the involutive algebra that facilitates the present analysis. 

However, in contrast to the analysis in [y, 0] that used the Moyal star product, which, as shown in 
0, stems from a particular choice of basis in the quantum Hilbert space, our analysis is carried out 
in a basis independent operator approach that bypasses the use of any star product (as explained in 
@) and any ambiguities that may result therefrom. Indeed, the analysis carried out here becomes 
problematic for the Moyal star product as the necessary positivity condition is not satisfied 9]. On 
the other hand, the power of the operator approach was demonstrated by solving for the spectrum 
of a spherical well potential in a non-commutative plane [Io| and the Coulomb problem in 3D non- 
commutative space [12j (with the coordinates satisfying a SU(2) algebra ). The current approach 
and that of @, 0] also differ on a more fundamental level in that the analysis in [(| 0] was essentially 
carried out in the classical configuration space T-L c (see section [TTJ) . Here we generalize the notion of the 
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Conncs' distance to the true physical states of a non-commutative system that are density matrices 
on the quantum Hilbert space, Hq, which has a natural tensorial structure H q = T-L c ® Ji* (see section 
HT| . This offers rather tantalizing possibilities in terms of the modification of the implied geometry and 
the possible statistical underpinnings thereof. 

The paper is organized as follows: In sectionUHwe briefly review the formulation of non-commutative 
quantum mechanics on a two dimensional non-commutative (Moyal) plane and in section lTOl we provide 
a brief review of earlier notions of geometry in quantum mechanics, which used the inner product on 
the Hilbert space to induce a metric on the manifold of quantum states. In section Hvl we show how the 
results of @ can be recovered from a spectral triplet defined on the classical configuration space and 
establish a link with the earlier notions of geometry discussed in section IIII1 In section [V] we proceed 
to introduce a spectral triplet on the quantum Hilbert space and introduce distances between mixed 
states. We also explore the possible statistical underpinnings of the emerging geometry. Section I VII 
summarizes and concludes the paper. 



II. QUANTUM MECHANICS ON THE NON-COMMUTATIVE PLANE 

In two dimensions the non-commutative plane is defined through the commutation relations 

[fi, £j] = itijd. (2) 

where e is the anti-symmetric tensor with t\2 = ^£21 = 1< We use a hat to emphasize the operator 
nature of these coordinates. One constructs standard creation and annihilation operators and b: 

b= x 1+ ix 2 ]= x x -i£^ 

and view the non-commutative plane as a boson Fock space spanned by the eigenstate \n) of the radial 
operator b%. We refer to it as the classical configuration space W c : 

n c = span{\n) = -L(tf) n \0)}. (4) 



This space plays the same role as the classical configuration space 1Z 2 in commutative quantum me- 
chanics. Next we introduce the quantum Hilbert space in which the physical states of the system 
and the non-commutative Heisenberg algebra are to be represented. This is taken to be the set of 
Hilbert-Schmidt operators over "H c and we refer to it as the quantum Hilbert space, Hg, 

U q = {ib : tr c (V>ty) < (5) 

where the subscript c refers to tracing over T-L c . Elements of "Hq are denoted by a round bracket l^) 
and the inner product is defined as 

(</#) = tr c (0ty). (6) 

We reserve f to denote hermitian conjugation on the classical Hilbert space, while \ denotes hermitian 
conjugation on the quantum Hilbert space. Note that H. q has a natural tensor product structure in 
that it can be viewed as the tensor product of H c and its dual, i.e., H q — T-L c ® H*. We can therefore 
also write the elements of T-L q in the form \ip,(f>) = We shall use this notation quite often 

below and refer to ip as the left hand and <f> as the right hand sector. 

In units of h — 1 a unitary representation of the non-commutative Heisenberg algebra (we use a hat 
to emphasize the operator nature of the coordinates and momenta and to distinguish them from their 
classical counterparts) 

[£i,Xj] = ieijO, (7) 



(8) 
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\PuPj}=0 (9) 

is obtained by the following action: 

Xi\i/>) = \xnl>), (10) 

Pi\1>) = \\ea[*iM)- (11) 

Here we use capital letters to distinguish operators acting on the quantum Hilbert space from those 
acting on the classical Hilbert space. For notational simplicity we also drop the hat notation as the 
capital letters distinguish these operators from their classical counterparts. It is also useful to introduce 
the following quantum operators 



P = P 1 + iP 2 , 

P t = P 1 -iP 2 . (12) 

These operators act as follow 

B\1>) = H), 

p N0 = -*yfl[M), 

P X W) = i\Jl\[tfM)- (13) 

The interpretation of this quantum system now proceeds as for a standard one. The only mod- 
ification required is that position measurement must now be interpreted in the context of a weak 
measurement (Positive Operator Valued Measure) rather than a strong (Projective Valued Measure- 
ment). The essence of the construction is based on the minimal uncertainty states on non-commutative 
configuration space, which are the normalized coherent states 



|0>, (14) 



where z = — }= {x\ + ix 2 ) is a dimensionless complex number. These states provide an overcomplete 
basis on the non-commutative configuration space. Corresponding to these states one constructs a 
state in quantum Hilbert space as follows: Introduce the non-orthogonal projection operators \z)(z\ 
on H c and define 

\z,z) = \\z){z\), B\z,z) = z\z,z), (15) 

which leads to the natural interpretation of {xi,x 2 ) as the dimensionful position coordinates. These 
states provide an overcomplete set on the quantum Hilbert space in the form [ill . Il3j 

l q = [ ™*\z,z)M{z,z\, = f) I (16) 

<— — ► 
where dg denotes differentiation to the left and d z differentiation to the right. The operators 

n z = ^\z,z)e^(z,z\ (17) 
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form a set of complete, positive, but non-orthogonal and unormalized projection operators [13[ , i.e., 

J dzdz n z = l q (<0|7T Z |V>) > 0,Vip, tt z t: w ^ S(z - w) , tt 2 z cx tt z . (18) 

They therefore provide a Positive Operator Valued Measure (POVM) that can be used to give a 
consistent probability interpretation by assigning the probability of finding the outcome of a position 
measurement to be (21,2:2), given that the system is described by the density matrix p, to be 

P{xi , x 2 ) = tr q (w z p) , (19) 

where tr q denotes the trace over quantum Hilbert space. In particular, for a pure state density matrix 
p = \if})(if>\, this probability is given by 

P(ii,i 2 )=tr q (7r,p) = (VMVO = 2^M*)/'^C#), (20) 

which differs from the commutative probability only by the presence of the Voros product fllj ] . 
In [l3|, [l4[ it was shown that these operators can also be written as follows: 

^z = -^z^2\z,n)(z,n\, (21) 

n 

where \n) is an arbitrary orthonormal basis in the classical configuration space, implying that this 
probability should be interpreted as the "total" probability of finding (21,22) as the outcome of a 
position measurement, which is insensitive to the value of the additional quantum number n labeling 
the right hand sector. 



III. INDUCED METRIC FROM THE HILBERT SPACE INNER PRODUCT 



It was realized some time ago that the inner-product on a generic Hilbert space % quite naturally 
induces a metric on manifolds of quantum states (unit rays In this section we briefly review 

this, however, instead of working with the quantum states as unit rays, which invariably restricts the 
formalism to pure states, we use the density matrix as the primary object that represents the states 
of a quantum system. This has the advantage of generalizing this notion of geometry to mixed states 
and it also ensures that the results are manifestly gauge invariant. This will also establish a link to 
Connes' notion of spectral distance in the next section. 

Let us therefore consider a n-parameter family of density matrices p(s) with s = (si, ....,s n ) € 1Z n . 
The density matrices have the usual properties of hermiticity {p* — p), unit norm (trp(s) = 1), semi- 
positivity and we assume that the map s 1— > p(s) is at least C . Realizing that the density matrices 
can be equipped with the natural inner product and trace norm on Hilbert-Schmidt operators, 

(A,B) = ti(A^B), \\A\\ 2 =(A,A), (22) 

the distance between two neighbouring states (density matrices) with coordinates s and (s + ds) can 
be defined as 

d 2 {l s + ds) = \\\5p || 2 = ^(d iP , d J p)ds l d^. (23) 

In the case of a family of pure state density matrices p(s) = \ip(s))(ip(s)\, \ip{s)) e "H it is a matter 
of straight forward verification that the metric tensor 

gij = -(dip,djp), (24) 

when re-written in terms of ip(s) and its derivative, yields [TBI. fl(| 



= Reidiifj, djip) - AiAj. 



(25) 
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Hero A- t = —i(rp,diip) represents the pull-back of the connection 1-form of the U(l) principal bundle 
onto the parameter space, and the exterior derivative of this connection 1-form A = Aids 1 gives the 
symplectic 2-form: a = dA = Im(diip,djifj)ds l A aV. This can be further expressed in terms of the 
covariant derivative Di = (di — iAA as gij = (Ditp, Djijj). 

Finally, d 2 (s, s + ds) can indeed be shown to be the smallest distance in the space of rays in the 
projective Hilbert space [l5[ 

d 2 (s,s + ds) = D 2 (R l p 1 ,R^ 1 ) = inf || V^e lQl - ^e"* 2 || 2 (26) 

ai,c*2 

where R^, represents the ray associated with the state ip. 



IV. SPECTRAL TRIPLETS ON CLASSICAL CONFIGURATION SPACE 

We start our analysis of the geometrical content of the construction in section |TT] by identifying the 
spectral triplet {A, W , D) where A = H q , % = "H c <8> C 2 , the action of elements of a G A on H is defined 
through the representation ir(a)(\ip), |</>)) = (a\tp) , a\cf>}) and the Dirac operator is defined as 



D = -i\ -A I n ■ (27) 




Here we have chosen the Dirac operator as in [6j to enable direct comparison with their results. 

States lo are positive linear functionals of norm one over A. Pure states play a rather fundamental 
role and is defined as those functionals that cannot be written as a convex linear combination of two 
other functionals. The Connes' spectral distance between two states is then defined by 

d(cu, lu') — sup |w(a) — ui'(a)\, 

B = {aeA:\\[D,7r(a)}\\ op < 1}, 

P|| op = supM. (28) 

1 101 1 

This is a rather difficult quantity to compute and therefore we first cast it into a more tractable form 
for the purpose of explicit computation and interpretation. For this we shall restrict our analysis to 
states satisfying the following conditions: 

• The states u, uj' are normal states (see ref [I?} for the definition). 

• The states uj and uj' are separately bounded on B, i.e., uj(a) < oo and uj' {a) < oo, Va G B. 



Let Vq = {a e A : \ \[D, 7r(a)]|| op = 0}, then the states uj, uj' are such that uj(a) — uj'(a) = , Va € 
V . 



The first two conditions are actually quite mild and essentially only implies that we restrict our 
analysis to states that can be represented by density matrices. Generally, the third condition places 
a restriction on the states for which our result for the Connes' distance holds, but under some rather 
generic conditions on the Dirac operator, that implies a certain irreducibility condition, it turns out 
that the result holds for all states satisfying the first two conditions. To show this we note that their 
is a close link between normal states and density matrices. Indeed, any normal state over A can be 
uniquely written as [l7[ 

ui(a) = tr c (p u a) (29) 

where p u is a hermitian, semi-positive operator on T-l c of trace one. Note that since tr c ( j o^) < tr c (/9 w ) = 
1, p u is a Hilbert-Schmidt operator and therefore p u G A = H q . If lu is a pure state, so is i.e. it is 
of the form p^ = for some \tp) G H c - This result can actually be understood quite easily from 
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the Hilbert space structure of A = H q as it follows simply from the Riesz-Frechet theorem, where the 
properties of p w derive from the properties of lu. 

Using this, and noting that our assumption u(a) < oo and ui'(a) < oo, Va G B implies tr c (p w a) < oo 
and tr c (p lAj ra) < oo, Va G B and setting dp = p — p 1 , we can write the third condition as 

• The states are such that 

tr c ((a»a) = [dp, a) = 0,Va G V . (30) 

We use the notation dp to indicate that we are interested in infinitesimal changes. Furthermore note 
that since tr c (dp 2 ) < 2(1 — tr c (pp')) < 2 from the semi-positiveness of p and p' , dp is also a Hilbert- 
Schmidt operator, as one would also have expected from the Hilbert space structure of these operators. 
Let us consider a Dirac operator of the form (|27[) 

D — *(co)' (31 > 

with the requirement that c, forms an irreducible pair on T~L Cl i.e., [c, a] = [c^,a] = =>■ a oc l c . It 
is easily checked that 7r(a)]|| op = => ||[c,a]|| op = , a]\\ op = =>■ [c,a] = [c',a] = => a oc l c . 
Since the identity element is not Hilbert-Schmidt, the set Vq contains only the null element, and 
the condition tr c ((dp)a) = holds for all states satisfying the first two conditions. For the choice 
c = b in (|27[) it immediately follows that the Dirac operator satisfies this condition since T-L c carries an 
irreducible representation of the oscillator algebra. 

Returning to the Connes' distance function, we can write 

d(p, p') = sup |tr c ((dp)a)| = sup |(dp,a)| (32) 

We are now in a position to determine the supremium on the right of Q32p . Let us introduce the 
orthogonal complement, Vq, of Vq. Note from (|3"U)) that Vq is a subspace of the space of elements 
orthogonal to dp. We can then decompose any a € A uniquely into its component a G Vq and its 
component eti G V (note that [D,ir(a)] ^ 0,Va G V ) . From (f3T)|) we can then write the Connes' 
distance function as 

d(p,p') = sup |tr c ((dp)a)| = sup \(dp,a)\, (33) 

aeB' aGB' 

B'={aeVj-:\\[D, ir (a)]\\ up <l}. 

The first step is to find a lower bound. To do this we note that the family of operators A = {a G 
A : a = Adp, < A < n-rg ^A^rn — } C -B'. Taking the extremal element a — rrrg ^t^tm — yields the 
lower bound 

«W)> nJ'tl'L ■ ( 34 ) 

|][-D,7TW)]||op 

Next we show that this is also an upper bound. Consider any a G B' and write it as a = ||a|| tr a 
where ||a|| tr = 1. Then we have from (|33| 

IH| tr ||[A^a)]||o P <i. (35) 

Since \\[D, 7r(a)]|| p > 0, this places an upper bound ||a|| tr < ^ with s = inf^gs' ||[Z?,7r(d)]|| op . 
We decompose a = cos(8)dp + s'm(6)dp ± (d) with dp and dp ± (d) unit vectors parallel and orthog- 
onal to dp. Note that both components are still in and that the latter depends on d. Since 
||[Z>,7r(a)]||op < |cos(0)|||[D,7r(dp)]|| op + | sin(0)|||[£>, 7r(o> ± (a))]|| op , with the equality holding for 
9 = 0,f we have s < s' = infagB'fl cos(6)\\\[D,ir(dp)]\\ op + \ sm(9)\ \\[D, ir(dp x (a))]\ | op ). Varying 
only the component of d along dp, i.e., varying 9, while holding dp_|_(d) fixed (geometrically this 
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corresponds to rotating a in the plane formed by a and dp), it is simple to see that the minimum is de- 
termined by the extremal points 9 = 0,^, depending on whether |[D,7r(d/3)]|| p or ||[£),7r(<ip_ L (a))]||op 
is the minimum (we can restrict ourselves without loss of generality to 9 6 [0, £]. Thus we conclude 

s < s' = min{||[£>,7r(dp)]|| op ,infaeB''(||[-D,7r(a)]||op)} = min{si, s 2 } with B" = {a e B' : (dp, a) = 0}. 
However, since for 9 = 0, § the equality holds in the triangular inequality, we must have s — s' . Since 
the right hand side of (f3"3"f vanishes for a € B", wc conclude that regardless of whether si or S2 is the 
minimum, 

\(dp,a)\ < \\dp\\ tl .\\a\\ tI < = mn lld fl m (36) 

si ||[D,7r(dp)J||op 

This establishes this value as an upper bound and therefore the supremium must satisfy 



d ^P')^ urn TTWi ■ ( 3? ) 

\\l D ,ndp)}\\ op 

Inequalities ([M]) and (I3T1) together yield the desired result 

This is a much more tractable form for Connes' spectral distance, which clearly has a strong resem- 
blance with (|23|) in the previous section, although the computation of the operator norm can still be 
tedious. In fact, in the examples we compute below it turns out that if we replace the operator norm 
with the trace norm it only modifies the distance function by a trivial numerical constant. Although 
a general proof that this is always the case lacks, it does suggest that defining a distance function, 
closely related to the Connes' distance function, by replacing the operator norm with the trace norm 
in psp is sensible and leads to simplification since the trace norm is often easier to compute. 

It is worthwhile highlighting the role of (|30|) in deriving (l38|) . If this condition does not hold, we can 
find a ao € Vq such that (dp, do) ^ 0. Since ([28| now places no constraint on ||ao||tr> it is clear that 
no upper bound can be derived for (dp, an) and it is in principle unbounded. 

Let us apply this to the example in Q where the distance between harmonic oscillator states is 
calculated. In the harmonic oscillator basis we can write a — J2 m n a m,n\m) (n\. The states are 
taken as Lu rn (a) = a m ,m with corresponding pure state density matrix p m — \m)(m\. We then take 
dp = p m +i — p m - The calculation is now straightforward: 



[D,n(dp)\ 




[b\dp] 
lb, dp] 



A simple calculation yields 

[b, dp] 1 * [b,dp] = 4(to + l)|m+ l){m + 1| + m\m)(m\ + (m + 2)|m + 2)(m + 2|, 

[6 f ,dp] f [b\dp] =(m + 2)|m + l){m+ 1| + 4(m + l)\m)(m\ + m\m -l){m- 1|. (40) 

As both of these operators are already diagonal, the operator norm, which is nothing but the largest 
eigenvalue can be read off exactly. This yields 



||[ jD)7r (^)]|| op = 2 A /H^tl) ,41, 



and 



^ + 1 ' W ) = V^TT)' (42) 
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which agrees precisely with the result of Q . 

Next we calculate the distance between two coherent states. We define the usual harmonic oscillator 
coherent states in H c given by \z >= e~^e zb |0) and consider the corresponding pure state density 
matrix p — \z)(z\. Then 

dp = \z + dz) (z + dz\ — \z)(z\ 

= (tf-z) \z)(z\dz+\z)(z\{b- z)dz + 0(dz 2 ,dz 2 ). (43) 

This can be written more simply as 

dp = dz\0)(l\+dz\l)(0\ (44) 

where we have introduced new bosons b = b — z, w = — z and the vacuum |0) = \z) for which 
6|6) =0, |1) = &t|o>. From this 

[b, dp] = dz\0)(6\ - V2dz\b)(2\~dz\l)(l\ (45) 

and 

[6,dp] f [b,dp] = dzdz\0}(6\ - V2(dz) 2 |0)(2| - V2(dz) 2 |2)(0| + 2dzdz\2}(2\ + dzdz\l)(l\, 
[b\dp] ] [b\dp\ =Mzdz\b){Q\ + dzdz\l){l\. (46) 
This yields, on diagonalising the first one, 



||[A7r(dp)]|U = (47) 



and 



20 

d(z + dz, z) = \ — V dzdz, (48) 
V o 

yielding the Euclidean distance. One can quickly check from (|4T)]) and (H6]l that replacing the operator 
norm in (|42|) and (|48|) by the trace norm one recovers, up to numerical factors, the same results. In 
what follows we shall, for computational simplicity, therefore often replace the operator norm with the 
trace norm and rather compute the closely related distance function 

^P')= W » c{ f m (49) 
||[D,7r(dp)]||tr 

In the examples above this distance function and the Connes' distance function only differ by a nu- 
merical constant. Whether this holds generally is an open question that needs further investigation. 



V. SPECTRAL TRIPLETS ON QUANTUM HILBERT SPACE 



In the construction of the previous section there is a one-to-one correspondence between states 
(points) \ip) in the classical configuration space and pure state density matrices p — acting 
on classical configuration space. These density matrices are also elements of the quantum Hilbert 
space corresponding to diagonal states \ip,ip) = ||V')(V'I)- However, there are many more physical 
pure states in the quantum Hilbert space of off-diagonal form \ip,<fr) = ||^)(^|). These states are in a 
one-to-one relation with the pure state density matrices p q = \i/),<f>)(ip,<f>\ acting on quantum Hilbert 
space and in a many-to-one relation with states (points) in classical configuration space, while only 
the diagonal density matrices p q — \?p, tf>)(tp, tp\ are m a one-to-one correspondence with states in the 
classical configuration space. Note that this is quite different to commutative quantum mechanics 
where there is a one-to-one relation between pure states density matrices and points in configuration 
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space. Specifically, in commutative quantum mechanics a particle localized at a point x is described 
by the pure state density matrix p = \x)(x\. On the other hand, as discussed in section [TT1 and further 
below, in non-commutative quantum mechanics a particle maximally localized at a point {x\,X2), 
without any prior knowledge of the right hand sector, is described by the mixed state density matrix 
p q = J2 n \ z ' n )( z > n \ ( z = - T i +ix2)- When the right hand sector is also specified the particle is described 
by a pure state density matrix p = \z, n)[z, n|, but note that there are infinitely many choices for the 
right hand sector given a point (£1,0:2) in classical configuration space. 

We would therefore also like to calculate distances between off-diagonal pure states and mixed states 
on quantum Hubert space. To do this we introduce a further spectral triplet (A,T-L,D) on quantum 
Hilbert space where A are the Hilbert-Schmidt operators on W q , T~L = H q ®C 2 , the action of elements of 
a G A on H is defined through the representation ir(a)(\ip), \(j>)) — (a\tp) , a\<p)) and the Dirac operator 
is defined as 




D = ~ i \h b n • (5°) 



Note that the Dirac operator defined here acts only on the left sector as it does not involve any 
momentum operators. 

Clearly, under the same conditions as in section HVl on the states, the analysis that led up to 
holds here with the replacement H c H q , i.e., 



trq(rfpg) 2 
||[D,7r(dp g )]|| p 



1/ f\ "qW«; /r-|N 



However, since the Dirac operator (JSOJ) only acts on the left sector, the condition (f5U|) becomes non- 
trivial. Indeed it is simple to verify that the general operators that commute with the Dirac operator 
are of the form 

T = ^2J2 a k,e\P,k)(p,i\ (52) 

p k,£ 

with \p) an orthonormal basis in classical configuration space and <Xk,i arbitrary complex numbers 
independent from p. In this case these operators span Vq. Generically (dp q ,T) ^ 0, (j30f is violated 
and we cannot expect (f5Tj) to yield the Connes' distance between all states. At best it will yield a 
lower bound, while the true Connes' distance may even be unbounded. However, (f5Tj) still applies for 
states such that (dp q ,T) = 0, VT and can be used to compute the Connes' distance function between 
such states. Below we apply ([5T]) in this context. Indeed, we shall compute (l5T|) regardless of whether 
the condition (dp q ,T) = is met or not and simply keep in mind that this will only yield the true 
Connes' distance when this condition holds. This is particularly useful in the context of a variational 
calculation as performed below. 

For an infinitesimal change dp q in the density matrix p q , a calculation paralleling the calculation 
leading to equation yields 

rn (a Mtrn (a w 2 ( [B^prf [B,dp q ] \ 

To start we compute the distance between the pure states corresponding to the density matrices 
p q (m, 4>) = |m, 4>)(m, <p\ with |m) a harmonic oscillator state and \(f>) an arbitrary state in the classical 
configuration space. Introducing dp q = \m + l,cj))(m + l,t/>\ — \m, (j))(m, (f>\, the calculation proceeds 
exactly as in section HVl Indeed, the state \4>) is purely a spectator and the result obtained is exactly 
the same as in (jl2")l . The same is true for the pure states described by the density matrix p q (z, <fi) = 
\z , (f)){z , <j)\ where \z) is a coherent state in the classical configuration space. Computing the distance 
between these states yields again the result (|48l) . independent of \cf>). The independence of these results 
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from the choice of the right hand sector specified by \(j>) is of course a consequence of the choice of the 
Dirac operator, which only acts on the left hand sector. In these cases it is also simple to verify that 
(dp q ,T) = 0, VT and thus this will yield the true Connes' distance function. 

We can use the same Dirac operator to compute the more general spectral distance between the 
states p q (m,<p) — |to,0)(to, 0| and p q (m + 1,0') — \m + 1,0') (m + 1,0'| where 0' is different and 
orthogonal to 0, i.e., (0|0') = ^<j>,<j>'- Introducing dp q = \m + l,0')(m + 1,0'| — |m, 0)(m, 0|, we get 

[B,dp q f[B,dp q ] = (m+1) | m+ l,0')(m + 1,0' | +(m + 1)8 ^ \ m + 1, 0')(to + 1, | +m | m,0)(TO,0 | 
+ (m + 2) | rn + 2, 4>')(m + 2, 0' | +(m + 1)<5 00 , | m + 1, 0')(to + 1, 4>' I 
+ (to + 1) | to + 1, 0)(m + 1,0|, (54) 

[B^dp^iB^dpq] = (m + 2) | m + 1,0')(to + 1,0' | +(m + 1) | to,0)(to,0 | +(m + 1)^, | to,0)(to,0' | 
+ (m + 1)5^/ | m, 0')(m, 0' | +(m + 1) | to, 0')(to, 0' | +to | to — 1, 0)(to — 1, | . (55) 

(56) 

It is easily seen that for = 0' we recover the earlier result (1421) . but for 7^ 0' we get 



||[£> )7r (dp g )]|| op = W^±^ (57) 



2$ 

This eventually yields d(p(m + 1, 0'), o(m, 0)) = \/ . In this case, however, (d/3 Q ,r) 7^ 0, Vr so 

V to + 2 

that this does not correspond to the true Connes' distance. However, taking this quantity at face value, 
it does show that d(p(m + 1, 0), p(m, 0)) < d{p{m + 1, 0'), p(m, 0)) and demonstrates that, despite the 
fact that the Dirac operator only acts on the left hand sector, the distance between two states depends 
on the right hand sector and, indeed, that the distance increases when the right hand sectors are taken 
differently Here we have taken the right hand to be a pure state that depends on the left hand, i.e., 
it changes from point to point. One may contemplate an even more general situation by taking the 
right hand sector to be a statistical mixture that changes from point to point, i.e., we can consider the 
following density matrices 

p q (m) = y^p n (m)\m,n)(m,n\, ^p„(m) = 1,Vto, 

n n 

P q ( z ,z) = ^2p n {z,z)\z,n)(z,n\, ^2p n (z,z) = 1, V2, (58) 

n n 

where \m) and \z) are as above, while |n) is an arbitrary orthonormal basis in classical configuration 
space and p n a set of position dependent probabilities. To clarify the physical meaning of these states, 
we compute the average of the radial operator in p q {m), which simply yields tr q (B\BLp q {m)) — to. As 
the fluctuations vanish, these are states localized at a fixed radial distance to. Similarly the averages 
of (B$ + B)/2 and i(B$ — B)/2 in the state p q (z, z) yield the real and imaginary parts of z, while the 
their fluctuations saturate the minimal uncertainty relation. This implies that these states describe a 
system maximally localized at the point (2:1,2:2). Note that the operators ir z introduced in (|2"Tj) are 
density matrices of this type as they can be interpreted as the density matrix for a particle maximally 
localized at z, but with complete ignorance of the right hand sector, i.e., all probabilities are equal. 
Indeed, computing the probability of finding a particle at point z° = (x® + ix^/VW, given that it is 
prepared in the state described by the density matrix it z with z = {x\ + 1x2) /V20, yields a Gaussian 
P(2;i,2;2) oc e~5ff (Oi-zJ +(x 2 -x 2 ) )_ ^ j g a j g0 f or reason that their use as a POVM leads to a 
measurement of position that disregards any information of the right hand sector. 

In these cases it is easily verified that if the probabilities p n are position independent, i.e., they do 
not depend on to or z, (dp q ,T) = 0, VT. Thus in this case (IBTI) will yield the true Connes' distance 
between these states. 

The choice of probabilities in (|58j) is, of course, an open issue. In equilibrium statistical mechanics, 
the choice of probabilities is dictated by equilibrium considerations, i.e., one maximizes the entropy 
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subject to constraints on averages, e.g., the average energy that leads to a Boltzman distribution. 
Here it is not clear which criteria should be used to determine the probability distribution, but two 
obvious possibilities present themselves. The first is to choose the probabilities such that the path 
length between two points is minimized. The second is to fix the probabilities from a condition of local 
thermal equilibrium, i.e., the local average energy is fixed and the local entropy associated with the 
density matrices (|58[) is optimized. The rest of this section explores the interplay between geometry 
and statistics from these two perspectives. Ideally one would prefer that the probability distributions 
and geometries that emerge from these two perspectives coincide. However, it turns out below that 
this only happens in the T — > oo limit. 

To facilitate the computation below, we use the modified distance function (|49]l , adapted to quantum 
Hilbert space, rather than (J2J), i.e., 



d{Pq,p'q) 

Setting dp q — p q (m + 1) — p q {m) one obtains 



d(m + 1, to) = 



\[D,7T(d Pq )}\\ tr - 

E„ (j£(»" + l)+j#(m)) 



2 (( 2m + 3K(m + 1) + (2m + l)p 2 n {m) + 2(m + l)p n (m + l)p n (m)) 



(59) 



(60) 



Let us first consider the first perspective and ask which choice of probabilities actually minimizes 
the distance function. Starting with (|60j) . the distance between two points mi and m/ is given by 



d(mf,mi)= J(to + 1,to). (61) 

m—rrii 

One easily concludes that the probabilities that minimize the distance must satisfy the equation 



where 





= A, Wi 




/ Pn(mi) ^ 




f A(m0 \ 




, A = 




Pnim+i) 




A(toj+i) 


\ Pn(m f ) J 




V A(to/) J 



(62) 



Pn = 



A(m) are the Lagrange multipliers imposing the constraints that the probabilities sum to one and A 
is the tri-diagonal matrix 



A = 



/ a{rrii) b(rrii) ^ 
b(nii) a(m i+ i) b(m i+ i) 



\ 



b(m f ) a(m f ) J 

The matrix elements of the matrix A are given by 

a(m) = V9(g( to) + g(m — 1) + (2m + l)(/(m) + f(m — 1))) , 
b(m) = -Ve(m + l)f{ TO + 1 ) . 

1 

g(m) 



VE n (( 2m + 3)i£(m + 1) + (2m + l)p 2 n {m) + 2(m + l)p n (m + l)p n {m)) 

/(m) = 1 En^t^+lj+PnW) 

2 (E„ ((2m + Z)pI (m + 1) + (2m + \)p\ (m) + 2(m + l)p n (m + l)p„(m))) 3/2 ' 



(63) 
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It is important to note that although A is a function of all the probabilities, it is independent of n. 
Using this we can sum over n in (1621) , which shows that A = Ap, where p is the column vector with 
equal entries I/O, being a cut-off for the sum over n. Substituting this back in (1621) shows that the 
probabilities are independent from m and n. In this case the distance function reduces to 

d(m + l,m)= (64) 
V60 v fn + 1 



which only differs from (l42l by a probability dependent global scale factor. In this case (dp q , T) = 0, VT 
and (|64| should be closely related to the Connes' distance, the only difference residing in the use of 
the trace norm in (|59p. rather than the operator norm. We expect this to only yield numerical factors. 

Next we explore the geometry emerging from the second perspective. To do this we introduce a 
local entropy associated with the density matrices (1581) in the usual manner 

S (m) = - J~]p n (m) logp„(m), 

n 

S(z, Z) = - ^Pn(z, z) \0gp n (z, z). (65) 
n 

A simple calculation shows that maximizing this local entropy leads, as above, to equally distributed 
probabilities. This, of course, corresponds to the T — > oo limit where only entropy plays a role and all 
states are assigned equal probabilities. 

If we associate an energy scale e n with the right hand states labeled by n and require the local average 
energy to be fixed on E(m), maximization of the local entropy yields a local Boltzman distribution 

-/3(m)e„ 

Z(fi{m)) ^ 

where /3(m) is the local inverse temperature. From this the distance function (|60|) can be computed. 
In particular, if we assume the local average energy and hence the temperature to be independent of 
m, we obtain 

- Ve ^zW) 1 

d(w + 1 ' m) = 7! W) v m + v (67) 

which again only differs from (|42p by a temperature dependent global scale factor. The same remarks 
as for (IMf regarding the relation to the Connes' distance are applicable here. 

This implied modified geometry for mixed states becomes even more explicit in the case of the 
continuous coherent state basis, which we consider next. 

The distance between the coherent state density matrices simplifies considerably with a slight redef- 
inition of the Dirac operator. We take the Dirac operator to be 

D ~ l ye{[B,27T0n z ] J (68) 

^ VnEnio,n)(i,ni o ( ™> 

Here we have again introduced the bosons 6 = 6 — z, ¥ = — z and the vacuum |0) = \z). Note 
that the Dirac operator is still insensitive to the right hand sector. Apart from the simplification that 
occurs with this choice of the Dirac operator, it is also a more natural choice in that it is based on the 
local vacuum |0) = \z) and local excitations W = W — z around this vacuum. The same remarks as 
made for the Dirac operator (|50l) apply to this Dirac operator. The general operators that commute 
with it are of the form 

f = X!I>M|P»*)(M (70) 

p k,l 
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As before, it is also straightforward to verify that for the density matrices (|58[). (dp q , f ) = 0, Vf if the 
probabilities are position independent. 

The computation now follows the same route as before and one obtains 



d(z + dz, z) = 1 



dp n (z) dp n (z)\ , ( f dPn{z)\ 2 ( dp n (z) s 2 



(71) 

This clearly exhibits a modified geometry. Assuming that p n is just a function of the dimensionful 



radial coordinate r = \j ^ one obtains 



d(z + dz,z) = { 2 p^ 2 +KM 2 ) V^TTw, (72) 

which is related to the Euclidean distance by a probability dependent conformal transformation. Note 
that since the probabilities are position dependent, this will not correspond to the true Connes' dis- 
tance. We can of course again compute the probabilities from the condition of local equilibrium. In 
the T — ¥ oo limit, the probabilities are obtained by maximizing (j65j), which yields equal, r-independent 
probabilities and subsequently, up to a global scale, an Euclidean geometry, which should be closely 
related to Connes' distance function. If we constrain the local average energy to E(r), the temperature 
becomes position dependent and the probabilities again follow a Boltzman distribution 

Pn{r) = wwr (73) 

yielding a conformal scale factor. If the local average energy is position independent, so is the temper- 
ature, and we obtain an Euclidean geometry with a global, temperature dependent, scale factor 



d(z + dz,z) = y/dr^+r^. (74) 



As the prefactor is bounded by < ^2Z{^) — \ > reacnm g hs minimum at T oo and its maximum at 
T — > this suggest, quite remarkably, that distances expand as the global temperature is lowered. We 
expect this distance function to be very closely related to the true Connes' distance as the probabilities 
are position independent and therefore (dp q ,T) = 0, VT. The only difference from the true Connes' 
distance is therefore the use of the trace, rather than operator norm, which we expect to yield only 
numerical pref actors. 



VI. CONCLUSIONS 



The formulation of non-commutative quantum mechanics on the Hilbert space of Hilbert-Schmidt 
operators naturally encodes the notion of spectral triplets and thus geometry, albeit in a sense stronger 
than in the conventional sense. Conditions under which a more tractable form of the Connes' distance 
formula apply have been derived. This form, and natural modifications of it, were used to make the 
geometry explicit on the level of the classical configuration space and the quantum Hilbert space. In 
the former case, as in commutative quantum mechanics, there is a one-to-one correspondence between 
pure state density matrices and points in the classical configuration space. In the latter case this 
no longer holds due to the tensor product structure of the quantum Hilbert space and the distance 
function depends on information encoded in the right hand sector, even though the Dirac operator 
only acts on the left hand sector. It then becomes natural to define a local entropy. The condition of 
local equilibrium then fully determines the distance function and, subsequently, the geometry. When 
the temperature is independent of position, Euclidean geometry results with a scale determined by 
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the global temperature. Remarkably the scale increases as the temperature is lowered, leading to 
expanding distances at lower temperature. 

There are a number of open issues worthwhile exploring further. This includes the consequence of 
the stronger technical condition on the involutive algebra imposed here and the issue of the modified 
distance function that invokes the trace rather than the operator norm. Another interesting point to 
explore further would be to investigate the class of Dirac operators for which (|5Tj) holds for all states 
and the implied associated geometries. 
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